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Abstract 

We study domain walls in two different extensions of super Yang-Mills characterized 
by the absence of a logarithmic term in their effective superpotential. The models, 
defined by the usual gaugino condensate and an extra field Y, give different patterns 
of domain walls despite both leading to the same effective limit for heavy Y, i.e. the 
Veneziano-Yankielowicz effective Lagrangian of super Yang-Mills. We explain the 
origin of those differences and also give a physical motivation for introducing the field 
Y. 



1 Introduction 



The dynamics of strongly coupled theories is a very interesting topic which has been the 
subject of detailed investigation for many years. It is now believed that the presence of 
supersymmetry (SUSY) helps in improving our understanding of these theories, by giving us 
a better insight into issues such as non-perturbative effects and the dynamics of confinement. 
In that respect SUSY would be not so much a phenomenologically relevant symmetry but a 
tool through which we could improve our knowledge of the structure of the theory. 

In particular supersymmetric Yang-Mills theory (SYM), the theory of gluons and their 
SUSY partners, the gluinos, offers a very attractive feature, namely the possibility of obtain- 
ing exact results despite it being a strongly coupled theory. In 1996 Dvali and Shifman [ij 
studied the formation of domain walls in this model and estimated the energy of such ob- 
jects when the walls were BPS-saturated (in which case also one half of the supersymmetry is 
preserved on the wall). Domain walls arise because SYM has T{G) degenerate vacua, where 
T(G) is the Dynkin index for the group G in the adjoint representation At some scale A 
the theory enters a strong coupling regime and gaugino condensates form, breaking the Z2N 
symmetry down to Z2. The gaugino condensate has been studied in detail for many years 
and by several groups using different methods 0. Given Aq, the gaugino field, we define the 
condensate as 

{XX) = {TiX^Xa) = A^e'^'' ,k = 0,...,N-l . (1) 

However if one wants to study SYM any further, it is necessary to find an effective description 
in terms of composite superfields. This was done in the early Eighties by Veneziano and 
Yankielowicz (VY) and the resulting Lagrangian is 
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where 



(few + h.c. 



(2) 

/C = {SSy/^ , (3) 



is the Kahler potential and 

>^ = §5 l^ln (^-^ ) -A^l , (4) 

is the superpotential, where S is the composite chiral superfield whose lowest component is 
the gaugino condensate 

S ^ ^.{Trico^co^) = ^(Tr(A.A")) + . . . (5) 

The structure of W is uniquely determined by the anomaly and the symmetries of the theory. 
Although this Lagrangian is perfectly adequate to reproduce the vacuum structure of SYM, 



^From now on we will consider G = SU(A^), therefore T(G) — N. 
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two major problems arise when one tries to use it to study the dynamics of the gaugino 
condensates. The first one is the lack of invariance of the superpotential which 
was pointed out by Kovner and Shifman a few years back0. The second problem is 
related to the presence of the logarithm in Eq. (|^), which results in the multivaluedness of 
the scalar potential. This problem was addressed (and solved) in the context of the linear 
multiplet formulation of gaugino condensation with a field dependent coupling p. In the 
chiral formulation, the concept of the glued •potential was proposed in Ref. [Q, and further 
developed in Ref. to cure both problems. However, it was pointed out in this latter 
reference how the presence of cusps (points where the effective potential is continuous but 
non-differentiable) between the different vacua of the model would constitute a problem 
when trying to construct domain walls. 

It therefore seems that the VY approach is somehow incomplete and we should be trying 
to incorporate new degrees of freedom in the effective description of SYM. Several proposals 
along this line have been put forward in previous years, and a certain degree of success 
has been achieved in applying them to the construction of domain walls. Promoting the 
scale parameter A to a dynamical field (the glueball order parameter) was one of them 0, 
which resulted in the construction of BPS-saturated domain walls in the large- limit [§ by 
exploiting the similarity of these models to a certain kind of Landau-Ginzburg construction. 
Another option is to incorporate chiral matter fields to the VY Lagrangian, something that 
was already considered in the Eighties by Taylor, Veneziano and Yankielowicz (TVY) |p. 
The idea is to study this model in its Higgs phase (i.e. for small values of the masses of 
the matter fields) and try to recover the VY limit when this mass parameter gets large. A 
lot of work has been devoted to following this approach in the context of building domain 
walls [|ry-|T^. In particular []TT| it is possible to construct BPS-saturated domain walls for 
any value of the mass parameter (and, therefore, obtain an effective VY limit) only when 
the number of matter flavours, A^/, is less than A^/2 and, moreover, the logarithmic branch 
of the superpotential is crossed if N j > 1. Although ways were found of circumventing 
these obstructions ||12|,0], the presence of the logarithm with its different possible branches 
is still not totally understood, and it would be desirable to find an alternative formulation 
where no logarithms were present other than in the limit were the VY theory was recovered. 
This is precisely the goal of this paper. A more ambitious objective would be to identify 
the extended models and new fields we introduce with already known physical systems. In 
particular, there has been a lot of work in recent years devoted to the connection between 



BPS-saturated domain walls and branes, starting with the work of Witten |T^, and we shall 
explore whether any connections are possible between our constructions, done in the low- 
energy field theory, and other studies performed in the context of string or M-theory ||16|| , 



and also in the presence of higher supersymmetries |17, 18 



In Section ^ we present two interpolating models, which we denote as Landau-Ginzburg- 
like, with their different interesting limits, one of which is precisely VY. In Section 3 we start 
by presenting a simple example, with A^ = 3, in order to illustrate the different aspects of the 
calculation, followed by a detailed discussion of more general results obtained in both models. 



^Note that, in terms of the S superfield, the invariance of the theory is Zj\j rather than the original 
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In Section 4 we motivate them from the physical point of view, and in Section 5 we study 
the large- expansion in the VY effective limit. We conclude in Section 6. We have also 
included a pair of Appendices, one explaining the numerical methods used to construct the 
domain walls, and the other supplying a detailed analysis of the BPS constraint equations. 



2 Interpolating models 

In this section we present two model supersymmetric field theories which interpolate between 
the Veneziano-Yankielowicz (VY) effective Lagrangian for supersymmetric Yang-Mills the- 
ory (SYM), and the so-called Landau-Ginzburg models. 

Let us consider the following superpotential which, from now on, will be referred to as 
Model I, 

MS,Y)=s(^^^ye-^/^^^-^ + ^^ , (6) 

with A and M two mass parameters. We will assume a canonical Kahler potential for the 
dimension one fields = S^^^ and Y, i.e. 

}C{S,S,Y,Y)^{SS)^/^ + YY . (7) 

Notice that the superpotential transforms as W — >■ e^^'^/'^W when S — >■ e^'^'^^^ S, which 
generates a symmetry in the Lagrangian. The supersymmetric vacua of this theory are 
given by solutions to the equations Ws — and Wy — 0, where 

Ws = (l-JV)(^)"e-''/«+''-i + I. , (8) 
There are N such solutions, 

= M{N - 1), S: = A^e^'™/^, (10) 

with a = 0, 1, . . . , — 1, which spontaneously break the Zn symmetry. In these minima, 
the superpotential is given by Wa = NS^. 

The mass eigenstates are most easily found in terms of the and Y fields. The mass 
matrix M is given by M = WW^, where W is the matrix of second derivatives of the 
superpotential 

w = ( y y ) . (11) 

In the vacuum (5°, Y"*), 
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In the limits of small and large A/M, one can straightforwardly show that 
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(13) 



(14) 



From now on we will work in M = 1 units. Then, for large values of A, y is a heavy 
field. We can therefore integrate it out by imposing Wy = 0, where Wy is given in Eq. (H). 
Substituting the resulting solution for Y in terms of 5* into the superpotential (^ we get the 
VY superpotential of Eq. (|), up to an unimportant constant factor, 



WyyiS) = S{\n 



S 



N 



A37V 



(15) 
= 



On the other hand, if we assume A <^ 1, then 5* is a heavy field, and we could impose Ws 
to obtain 

Wefr(r) = ATA^^ ( ^ j e-^/^+i-V^ . (16) 



The resulting potential has a single minimum at F = — 1, as is expected, which means 
that the trajectories described by Y start and end up at the same point. Moreover, it can 
be checked that, in the limit A — > 0, the field barely moves from that point. Therefore it is 
a good approximation to freeze Y at its vacuum value, to obtain an effective superpotential 



LGl 




+ N 



(17) 



The effective theory now has a form similar to a Z^r symmetric Landau-Ginzburg model 
(hence the notation), although it is not usual to see negative powers of the field. 

Alternatively, one can consider the following superpotential (our, from now on. Model II), 



WiS,Y) = S 



^ ^ Y/M-N-l _ x_ 

A3 j M 



which can be obtained from Eq. (^) by just changing (A^ — > 
are also A^ minima 

Y^^ = N + 1, = ASe^^-'^/^ 



-A^, Y —>■ —Y ). Hence, there 



(19) 



with a = 0, 1, . . . , A^ — 1, which spontaneously break the Zn symmetry. The superpotential 
in these minima is given by W* = ~NS^. It is easy to check that the same hierarchy between 
the S and Y masses appears in the large/small A limit. On top of that, there is also a Zn 
conserving minimum given by [Y^ = 0, S** = 0). Whereas for large A we get the same VY 
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limit of Eq. ([15|) , at small values of A the model approaches a conventional Landau-Ginzburg 
model 

>Vlg2 = 5(^(^) -iV-lj. (20) 

Hence we have come up with two different two-field models which, in the limit A ^ 1 
coincide with the VY Lagrangian, without the ambiguity of a multivalued logarithm and, 
as far as solutions to the field equations for S go, become Landau-Ginzburg models in the 
hmit A < 1. 



3 Construction of the domain walls 

Let us study the domain wall spectrum of these models. As we are considering static solu- 
tions, we wish to minimize the static surface energy functional 

^ = \j dz [JCssd.SdJ + Kryd^Yd^Y + V{S, Y, S, Y)) , (21) 

where JCjj = d^JC/dJdJ, J = S,Y is the corresponding Kahler metric. From now on we 
assume that the wall spreads along the xy plane and, therefore, the profile is calculated along 
the z direction. The scalar potential V is given by 

V = )C-l\Ws\^ + IC-\.\Wy\^ . (22) 
As is well known [D!^|2D|1, this can be rewritten as 



a 



where 7 is an arbitrary constant. Consider a domain wall interpolating between vacuum a 
and vacuum b. If we assume that it is BPS, it will be described by the following equations 

JCssd.S = e'^^ , 

(24) 

where now e"*'''"'' = {Wb — Wa)/|Wfe — Wa|, and Wa(b) is the value of the superpotential at 
vacuum a{b). One can easily show that, in that case, the energy functional saturates the 
BPS bound 

a=\Wt-Wa\ , (25) 

and 

9,(6^^""^) = V . (26) 
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Hence the solution to the BPS equations (if any) interpolates between Wa and Wb following 
a straight line in the W-space, obeying the constraint 

Im(e^^W(5,F)) = constant . (27) 

This constraint equation is the key to classifying the solutions to the BPS equations. 
From now on we shall adopt the following parametrization for the fields S and Y 

S = \S,\R{z)e'^^'^ , 

(28) 

Y = Re [Y]{z)+ilm [Y]{z) , 

where R and f3 are real functions. It is therefore easy to see that, when going from a vacuum 
a situated at z = — oo, to another one, b, which is k neighbours away (i.e. b = a + k), 
situated aX z = +00, the following boundary conditions apply 

• for the S field, 

R{-oo) = i?(+oo) = 1, Pi-00) = 0, p{+oo) = 2Txk/N 

• for the Y field. 

Re [Y]{-oo) = Re [Y]{+oo) = ^ 1, Im [Y]{-oo) = Im [Y]{+oo) = 

while the phase 7 is given by 7 = =F7r/2 — nk/N depending on which model (I or II) we are 
considering. We will take k in the range < < A^. Notice that there are two possible 
domain walls connecting two given vacua corresponding to the two paths defined by k and 
N-k. 

3.1 An example: N=3 

Before explaining how the different effective limits emerge, it is interesting to work out an 
exphcit example. Let us consider the simplest case, the extension of SU(3) super Yang-Mills. 
Since there are three vacua, we can always connect any two of them through a domain wall 
that interpolates between nearest neighbours. We want to find these domain walls and 
study their BPS character. In order to do that, we have to minimize the energy of a static 
configuration with the required boundary conditions. Our problem involves two mass scales 
and, therefore, the discretization has to be done in an optimal way to include them both. 
The details of the minimization algorithm and our convergence criteria are given in Appendix 
A. On the other hand, if we assume that the domain wall is BPS-saturated, we can try to 
solve the corresponding first order differential equations given by Eq. (p^). To do that, we 
have to find the appropriate (if any) initial conditions for the fields at the centre {z = 0) 
of the domain wall. Using a symmetric ansatz allows us to fix the phases of the fields at 
this point. We are then left with two parameters, the absolute values of the fields {R and 
Re [Y]) but, since they have to obey the constraint Eq. (P7|), only one of them is free. If 
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Figure 1: 

(a) Argand diagram of the profiles of the gaugino condensate S for TV = 3 and several values of A in Model I. 
The continuous line represents the effective VY limit; (b) Argand diagram of the profiles of the Y field for 
the same model and values of the parameters as before. 



there is a BPS domain wall, we can use the usual overshooting-undershooting method to fix 
that parameter, and then obtain the profile. 

We have used both methods. Of course, minimizing the energy gives more information 
since it allows us to explore also non-BPS domain walls. We have found that there is a BPS 
domain wall for all values of A, irrespective of which extension we are using. The domain 
wall path in S and Y space for the first model, defined by Eq. (|l3), can be seen in Figs la 
and lb respectively, where we plot the Argand diagrams for different values of A. Notice 
that they get stabilized both in the small and in the large A limits. 

As we will see, the same applies to the profiles when we plot them as a function of the 
relevant {z or Az ) variable. 

Let us first focus on the large A case. The profile for S is made out of two branches, 
z > and 2; < 0. According to our theoretical analysis of Section 2, this configuration 
should be also described by an effective VY Lagrangian. In fact, we have checked that it 
can be obtained by solving the BPS equations for the VY Lagrangian with an appropriate 
prescription for the logarithm. In particular, the z < {z > 0) branch is obtained using 
In 5*^ = ln(A'^i?)^ + A^i arg(5') (— 27ri). The profile for Y can be divided in three regions. 
The central one corresponds to the core, where 5* remains almost constant but Y changes 
very quickly, as can be seen in Figs 2. Notice that the width of this core goes to zero when 
we increase A. Moreover, the Y profile can be calculated from the 5* one. This is because, 
as expected in this limit, the condition Wy = is satisfied almost everywhere (it is only 
violated in the region where the phase is forced to change almost in an abrupt way to connect 
the two S branches, as can be seen in Fig. 2b). 

In the small A limit, the Y field is frozen to its value at the minimum, =p 1 depending 
on the model we are considering (as can be clearly seen in Fig. lb). We have checked that 
the profile for the S field coincides with the corresponding BPS domain walls in the effective 
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Figure 2: 

(a) Profiles of R (absolute value of S) and |y| as a function of Az for N = 3 and large values of A in Model I 
(solid: A = 100, dot-dashed: A = 20, dashed: A = 10); (b) profiles of Np/{2n) and arg(r) for the same 
model and values of the parameters as before. 

Landau-Ginzburg models of Eqs (p!7| , pOD . This is illustrated in Figs 3, where we show both 
the absolute value R (Fig. 3a) and the phase [3 (Fig. 3b) of the S field in Model I as a 
function of the rescaled position, Az, for small values of A, as well as its theoretical A ^ 
limit. 

We have therefore found a uniparametric family of BPS domain walls for the case (A^ = 
3, /c = 1). The theory describing these walls interpolates between the super Yang-Mills 
Veneziano-Yankielowicz model (with some given rules about how to glue the different sectors 
of the potential) and a Landau-Ginzburg model. On the other side, we have also analyzed 
the [N = 3,k = 2) case. There are no BPS domain walls with these boundary conditions (see 
Appendix Q), in other words the domain wall connects the two vacua through the shortest 
path in S space. 

Let us now try to extend this result for higher values of A^, where a richer pattern of 
domain walls is possible. Notice that BPS domain walls have to satisfy the constraint given 
by Eq. (|27|) . This equation involves two fields which, in the limits explained before (A — > oo 
and A — > 0), should split into two simple -more restrictive- equations, since the dynamics 
are controlled only by one field. As we will see, the compatibility of these equations is the 
key to understanding the pattern of domain walls for arbitrary values. 

3.2 General analysis: results 

When considering domain walls in models with A^ > 3, the possibility of having different 
types of walls, connecting the different vacua, arises. We first consider complex walls, which 
connect different vacua with 15*1 = A^, and briefly discuss real walls (connecting the vacuum 
in Model II at 1 5*1 = to a vacuum at 15*1 = A^) at the end of the section. 

For the moment we shall exclude the cases where k = N/2 (with A^ even). They corre- 
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Figure 3: 

(a) Profiles of R (absolute value of S), as a function of Az, for iV = 3 and small values of A (dotted: A = 2, 
dashed: A = 1, dot-dashed: A = 0.5, sohd: LG limit A ^ 0) in Model I; (b) profiles of 7V/3/(27r) for the 
same model and values of the parameters as before. The A = 0.5 curve is not shown explicitly as it is 
indistinguishable from the LG limit on the graph. 



spond to 5'(— oo) = — ^(oo) and give a constraint that includes 5 = as a possible point 
in their trajectory, where we have seen that the two models we are considering (I and II) 
have radically different behaviour. For a detailed study of these cases the reader is referred 
to Appendix 0. 

In our numerical searches for both BPS and non-BPS solutions, we have found just one 
domain wall solution interpolating between two given vacua. Our conclusions from these 
results and from our analysis of the BPS constraint equations detailed below are that, for 
finite A, 

• the domain walls found in Model I are BPS-saturated only for k = 1; 

• in Model II we have found BPS-saturated walls for any value of k in the range 1 < 
k < N/2. This actually covers all possible situations given that, for any A^, there is 
always a path with k < N/2 that links any two chirally asymmetric vacua. 

These statements are proved in Appendix |B] by going to the limit A — > and appealing to 
the continuity of the solutions as a function of A. Moreover they are illustrated in Figs 4 
for A^ = 5 and k = 2, where we plot, in Fig. 4a, the Argand diagram for the superpotential 
Eq. (|^) of Model I. As we can see the larger A is, the more the curve tends to a straight 
line, which corresponds to a BPS-saturated solution. Therefore, it looks likely that, in the 
A — >^ oo (VY) limit, these domain walls of Model I will turn into BPS ones. 

In Figs. 4b, c we turn to Model II and plot the profiles of the fields R (Fig. 4b) and Re [Y] 
(Fig. 4c) as a function of Az, for several values of A. The behaviour of the different curves 
is very similar to that of the A^ = 3, A; = 1 case plotted before, and we have also appended 
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in Fig. 4b the effective A (Landau-Ginzburg) and A ^ oo (VY) limits, which we shall 
discuss in a moment. 

Therefore we have seen that the pattern of domain walls for A; > 1 is totally different in 
Models I and II, although it looks like both are able to give BPS-saturated domain walls in 
the VY effective limit. In order to clarify this point further let us turn to the study of the 
constraint equation, Eq. (p^). 

Let us consider the walls connecting the vacua labelled by integers a and a + k (i.e. 
S{+oo) = S{—oo)e'''^'^^/^). Using the parametrization of Eq. (^), the constraint function 
for Model I is given by 



Ci{S,Y) 



j^l-N^-Rc [Y]+N-l 



COS 



[1-N)(3- 



Im [Y] 



+ RRe [Y] cos{P-^\ - Rim [Y] sin (p-^ 



A^ cos ( — 



(29) 



and, for Model II, 

Cn{S,Y) = i?^+ie^^[^l-^-^os 
-RRe [Y] cos ( /? 



7rk 
A^ 



+ Im [Y] 



nk 



Tck 



+ i? Im [r] sin ( /5 - — ) + A^ cos 



Trk 



(30) 



The full constraint equations Eqs (|29|j30|) are too complicated to study in general, but we 
can still gain some information about how the solutions interpolate between domain walls 
in the Landau-Ginzburg models and the VY limit by studying the corresponding constraint 
equations at 2; = 0. 

We are considering symmetric walls therefore, at the origin, = nk/N and Im [Yq] = 0, 
where we use the subscript to denote fields evaluated at z = 0. For Model I the constraint 
equation takes the form 



C,{Ro,Yo) = Rl-^e 



l-N -Yo+N-l, 



-IV + RqYo - N COS 



and, for Model II, 



Cu{Ro,Yo) = R^+'e 



N+lYo-N-l, 



RoYo + Ncos ( — 



Trk 



7ck 







(31) 



(32) 



In these equations Yq is real. 

We can now plot the relationship between Rq and Yq for both models, as a function of 
A^ and k. Falling on the constraint curve is necessary, but not sufficient, for the existence 
of a BPS solution. Figure 5 shows the curve given by Eq. ( |3TD for A^ = 5 and k = 1 (a) 
and k = 2 (b), and that given by by Eq. (|3^ for A^ = 5 and k = 1 (c) and k = 2 (d). 
The horizontal lines represent the value of Yq for which the corresponding Landau-Ginzburg 
(A — >■ 0) limit is reached {Yq = A^ — 1 in the first two plots and Yq = N + 1 in the last two). 
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Figure 4: 

(a) Argand diagram of W in Model I, N — 5, k = 2 for different values of A (solid: A ~ 100, dot-dashed: 
A — 20, dashed: A = 5); (b) profiles of R (absolute value of S) as a function of Az for = 5, A: = 2 and 
different values of A in Model II (dot-dashed: A = 5, dashed: A = 20, dotted: A = 40). The two limits are 
also plotted: VY A oo (thin solid) and LG A ^ (thick solid); (c) Profiles of Re [Y] as a function of Az 
for the same model and values of the parameters as before. 
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Figure 5: 

(a) Contour plot of the constraint equation at z = as a function of Rq and Yq (which is real) for = 5 
and fc = 1 in Model I; (b) same for k = 2; (c) same for Model II, fc = 1; (d) same for Model II, k ~ 2. The 
horizontal lines represent the corresponding values for Yq at the LG (A —>■ 0) limit, whereas the vertical ones 
represent the values of Rq in the VY (A ^ cxd) one. The annotated symbols indicate the BPS-saturated 
solutions found in each case. 
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Note that, in all but Fig. 5b, the constraint equation intersects the A — > horizontal line 
at either one (Figs 5a, c) or two (Fig. 5d) points. The intersection at larger Rq for Model II 
(A^ = 5, k = 2) does not in fact correspond to a BPS wall, as is shown in Appendix 

Let us consider now the VY (large A) limit, where the dynamics near the minima are 
governed by the light field, S. The Y field just follows the evolution of the 5* field, as dictated 
by the Wy = condition, F = A^ ^ 1 - ln{S^/A^^), with the minus (plus) sign taken for 
Model I (Model II). However, this condition must be violated in a region near 2; = 0, as it is 
not consistent with Im [Yq] = 0, /3o = nk/N. We denote the width of the region by A, such 
that (with a suitable choice of the branch of the logarithm) 

Y = NTl-\n^ z<-A , 

(33) 

Y = A^ ip 1 - In -^777 + i2'Kk z> A . 

As we will see, in that region R is approximately constant but Im [Y] + Nj3 changes very 
quickly to connect the two Wy = branches. We would expect the width of the region to 
be governed by the larger of the two masses in the theory, given in Eqs (0), in which case 
it vanishes as A~'^ in the limit A — 00. 

In this limit, and using Eqs (^), the constraint equation becomes, 

i?o(l-lni?o)=cos(^^^ . (34) 

This equation has two solutions, one with Rq < 1, and the other with i?o > 1- Our analysis 
in Appendix |B] shows that there is, at most, only one BPS wall for finite A, and that it 
must be the solution with smaller Rq. It follows that, in the A — > 00 limit, only one of 
these solutions to Eq. (|3^ ) corresponds to a domain wall, which is represented in Figs 5 by a 
vertical line, and denoted by Rqvy- Note also that the equation is identical to that obtained 
in Ref. [IT for the effective VY limit, with k replaced by Nf, the number of matter flavours 



introduced in that case to study the model in its Higgs phase. There the effective VY limit 
was obtained when m ^ 00, with m the mass of the matter condensates. Again, in that 
case, only the solutions with Nj/N < 1/2 (i.e. Rq < 1) would yield meaningful domain 
wall proflles, placing a restriction on the amount of matter one could consider. Here, as 
we mentioned before, the fact that it is only possible to deflne the effective VY limit for 
k < N/ 2 does not prevent us from building domain walls between any two neighbours (with 
the exception of A^ even, k = N/2). 

Once we have deflned the effective VY limit, let us go back to the discussion of our 
results, for which we focus again on Figs 5. It is noticeable how, in all but Fig. 5b, not only 
do the constraint curves intersect the Landau-Ginzburg (A — 0) limit, but also they connect 
continuously both the Landau-Ginzburg and VY (A — > 00) limits. These are precisely the 
curves along which we have found BPS-saturated solutions at flnite A, which are indicated 
by the annotated symbols, whereas the absence of such solutions is obvious in Fig. 5b, where 
the constraint curve does not connect the two limits. There is, therefore, a straightforward 
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criterion to decide whether a particular model would yield BPS-saturated solutions for finite 
A, and that is given by the fact that the corresponding constraint equation interpolates 
continuously between the two existing limits. In the particular case of Fig. 5d, where there 
are actually two branches of the constraint equation that link the VY and Landau-Ginzburg 
limits, it is shown in Appendix ^ that only one of them gives BPS domain wall solutions. It 
is further shown, in that same Appendix, that the shape of the constraint curves illustrated 
for N = 5 and /c = 1, 2 in Figs 5 are examples of generic classes distinguished hj k < N/2 
and k even or odd, and hence that our statements about the number of BPS solutions in the 
two models can be extended to all A^, k. 

Finally, as mentioned at the beginning of the section, there is the possibility, in Model II, 
of constructing so-called real domain walls, which interpolate between the chirally symmetric 
vacuum at {\S\ =0, Y = 0), and one of the A^ vacua with \S\ = A^, Y = N + 1. Although 
the chiral vacuum is certainly a supersymmetric minimum of the theory, in the sense that 
the scalar potential vanishes there, the Kahler potential (|^) gives a metric which is singular 
at that point, so the status of the chirally symmetric vacuum is ambiguous. Nevertheless we 
have verified numerically that solutions exist in the model for our choice of Kahler metric. 



4 Discussion, physical motivation, relation to previous 
work 

Once we have shown that both these two Landau-Ginzburg models admit domain walls, 
and also have a well-defined effective VY limit when A — oo, we can proceed to discuss 
what the possible physical meaning of the extra Y field should be. It was already pointed 
out years ago that the VY approach to SUSY QCD was neglecting degrees of freedom 
which would be as relevant as the gaugino condensate upon which the VY approach was 
developed. For example, spin zero glueballs which, in principle, should not be heavier than 
gluino condensates do not contribute to that effective action. Therefore it would be justified 
to attempt a different construction which would embed the VY model. This was done 
in Ref. by introducing a real tensor superfield U which contained gluino- gluino, gluino- 



gluon and gluon-gluon bound states, and a new U -dependent term to extend the VY effective 
action. Therefore this approach preserves the original logarithm and we do not believe that 
it could be directly compared to our Models I and II. 

As was already mentioned in the introduction, the problems associated with the presence 
of the logarithm in the original VY model were addressed in Ref. in the context of the 
linear multiplet formalism. In that case the authors were focusing on the study of gaugino 
condensation in the presence of a field dependent coupling, which brings into the theory a 
new scalar field whose vacuum expectation value induces the gauge coupling. By studying 
and proving the duality between the linear and chiral formulations of the problem, they 
were able to amend the VY chiral potential in order to eliminate ambiguities coming from 
the imaginary part of the gaugino condensate field {u in that case). Again, we have not 
been able to find a direct link between this approach and ours of extending the theory by 
introducing a second field Y but, nevertheless, we believe that exploring and constructing 
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domain walls in this now unambiguous formulation would be a very interesting exercise. 

On the other hand, our Y field certainly has the form of a field dependent coupling, as 
the effective action contains YS\f- If we accept that F is a coupling field, then it is tempting 
to interpret the extra piece, proportional to e^^, as the leading term of an instanton-induced 
contribution. 

Finally we can comment about the different applications of BPS domain walls in the 



context of string theory. According to 't Hooft's conjecture p2[, large- iV gauge theories 
should exhibit a phase described by perturbative strings. It is also believed that, in this 
limit, BPS-saturated domain walls are effectively branes on which the QCD string can end. 
Triggered by this, much work has been devoted to studying the connection between domain 
walls and branes in recent years. For example a construction of BPS domain walls was carried 
out in Ref. by exploiting the similarity of the super QCD Lagrangian (understanding by 
this the VY model amended to restore the Zjq symmetry 0]) with that of the A]^ Landau- 
Ginzburg model, in the large- limit and for walls connecting adjacent vacua. In fact, there 
it was even attempted to identify the field content of this super QCD Lagrangian with those 
fields appearing in the supersymmetric formulation of the membrane action. 

There has been also a lot of work on branes and strings at conifold singularities. Different 
dualities are supposed to be related to a transition in the geometry. This has been recently 



suggested by Vafa [jT6[ in the Chern-Simons/topological strings duality framework. He 
considers Type IIA strings in a non-compact Calabi-Yau (CY) 3-fold geometry that includes 
a conifold and adds A^ D6-branes wrapped over the in the complex deformed conifold. 
That generates a A/" = 1 supersymmetric SU(A^) theory. The relevant chiral superfields 
are 5", proportional to the one defined in Eq. (H), and a new chiral field, y, whose lowest 
component is given by the volume of the (real part) and by the vev of some 3-form 
(imaginary part) that will play the role of 9 angle. At lowest order, the superpotential is 
given by W = SY . The new field, Y, triggers the non-perturbative effects. Vafa proposes 
the following corrected superpotential 

W = \ SY + iN''e-^/'' , (35) 
A s 

where is the string coupling . If we integrate out the Y field, by using Wy = 0, we end up 
with the VY superpotential for the condensate. Notice that this superpotential shares two 
properties with our models: there is a linear, SY , term and an exponential term involving 
the new field, Y . The combination of these two facts gives the logarithmic dependence in the 
effective Lagrangian describing the gaugino bilinear. Therefore we believe that our models 
could shed some light on the role played by the fields that couple to the gaugino condensate. 
These fields are very important if we want to understand the core of the domain wall. 



5 The large- limit 

We have seen in the previous examples how to calculate BPS domain walls for the VY 
Lagrangian in the cases where they exist in the extended theory. Let us assume that the 
complete theory supports these BPS domain walls. Now we can study how they behave for 
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large values of N (see also Refs This is particularly important because, in that 

limit, these domain walls are thought to be the field theory realization of branes. 
We want to solve the BPS equation 

qN 

]CssdJ = e^''lnj^ , (36) 

for the {N,k) case (i.e. 7 = 7r/2 — nk/N), and large N values. In order to do that, it is 
useful to define the parameter e = Then the constraint for Rq, Eq. (|34D, becomes 

cose-i?o(l -lni?o) =0 . (37) 

Notice that the allowed values of k will depend on the details of the theory that generalizes 
the VY Lagrangian. For example, in Model I only k = 1 will be allowed^, whereas in Model II 
we can consider any value of k. In any case, assuming k <^ N, we can expand Ro in terms 
of e to obtain its large- limit 

Ro = l±e+'^ + Oie') . (38) 
6 

Then IS"! deviates from its value at the minima by 0{e). Since the phase of the condensate 
is also C(e), we will adopt the following parametrization 

S{z) = A%1 - tp{z)) e^^'^^") , (39) 

with p{z),ip{z) real functions of 0{1). Substituting the previous expression into the BPS 
equation, Eq. (p6D, we get 

p{z) = e^^/r + (e^^/r - -ei«^/r)e + 0(e2) , 

6 

(40) 

^{z) = e9^/^ + 0(e2) , 

for the z < branch (and the 'symmetric' one in the z > branch). In these equations 
r = ^ fixes the domain wall width. In order to determine the N dependence of the relevant 
quantities (i.e. width and tension of the wall) in the large-iV limit, it is standard to rescale 
the fields so that the Lagrangian scales like A^. For example, Eq. (|T|) will appear as 

(AA)fc = (TrA'^A,)^ = iVA^e^t'^ . (41) 

If we calculate now the domain wall tension when interpolating from vacuum a to vacuum 
a + A; we get 

<y = \Wa+k -Wa\= Ar2A32sine , (42) 

■^A similar result was found in Ref. j2^, where the large- limit is mapped into a Landau-Ginzburg 
model. 
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i.e. for k = 0{1) the tension is expected to scale as A^. Notice that there are two different 
ways to calculate the domain wall energy: using the two-field theory or the effective VY one. 
The validity of Eq. ( ^21) relies on the differentiability of the domain wall profile, so it can 
be only used in the two field theory. In the VY model, one should split the energy into the 
contributions coming from the z < and 2; > branches. In that case we get 

avY = N^A^2\ sin e - eRo\ . (43) 

This result coincides with the one obtained by using squark matter fields (see Ref. |T^) as 
the extra fields providing regular domain walls. This is due to the fact that, in both cases, 
the same prescription is used to calculate the imaginary part of the logarithm. We have seen 
how this prescription naturally emerges when integrating out the Y field (see Eq. (0)), as 
it was also shown in Ref. for the squark matter fields. Notice, however, that for A; > 1 



the squark fields, as described by the TVY Lagrangian, do not provide regular BPS domain 
walls for high enough values of their masses. This means that, in these cases, the extra fields 
that enter the logarithmic term in the effective potential can not be ignored. They play a 
relevant role and are excited in the sub core of the domain wall. 

On the other hand, the difference between a and ayy gives us information about the 
energy that is stored in the Y subcore, i.e. the spatial region in the large A limit where the 
phase of Y changes quickly. Since a ~ A^A^ and ayy ~ A'^, we conclude (as first shown in 
Ref. [14]) that, in the large- iV limit, most of the energy concentrates in this subcore. 



6 Conclusions 

We have studied the existence of BPS-saturated domain walls in two extensions of the 
Veneziano-Yankielowicz effective Lagrangian formulation in order to study SUSY Yang- 
Mills theories. These are motivated by the existence of several problems associated with 
the presence of a logarithm of the gaugino condensate field in the effective potential of the 
VY formulation, which causes ambiguities when studying dynamical questions such as the 
formation of domain walls. Our extensions, which we denote as of the Landau-Ginzburg 
type, avoid such problems by introducing an extra chiral superfield, Y, and a logarithm-free 
interaction which, in the limit of heavy Y, result in the effective VY model. Within this 
framework we have constructed domain walls which turn out to be BPS-saturated in one of 
the models but not in the other one. In any case, both seem to give BPS-saturated walls in 
the effective VY limit. 

We have also studied the limit of light Y field, in which the models proposed approach 
Landau-Ginzburg models with which one can construct well-defined BPS-saturated walls. It 
turns out that it is possible to define a criterion to decide whether our particular extensions 
of VY admit BPS-saturated walls at any point in between the two (VY and LG) limits: this 
is given by whether the BPS constraint, associated with the corresponding BPS equations, 
is able to interpolate continuously between both limits. Our results can be proven to hold 
for any values of (number of colours) and k (neighbours between which we interpolate). 

Finally we comment on the possible physical significance of the extra field Y, and the 
connection of our approach with other proposed models existing in the literature. Although 
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there is no direct evidence so far, the similarity of this model with some domain wall con- 
structions done in the context of string theory is quite remarkable, and it deserves further 
investigation. 
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A The Newton method for finding solutions 

The numerical solutions to the field equations were found with Newton's method, also called 
the Newton- Raphson method (NR). This was implemented by a Fortran 90 program. 

We wish to find extremal energy configurations for the discretized energy functional of 
the domain wall, E{fA), where /a denote the degrees of freedom: here, the values of the 
complex fields and their conjugates at each point on the lattice. 

The Newton-Raphson algorithm consists of iterating the update 

KdfAdfB J dfB 



In order to solve the matrix equation our program calls linear algebra routines zgbtrf 
and ZGBTRS from the Silicon Graphics implementation of the BLAS library. 

It is necessary to use a rescaled lattice because of the presence of widely differing mass 
scales when A differs from 0(1). The rescaling used takes the form 

s = -(tanh(mia;) + tanh(m2x)) , (45) 

where x is the physical position, s is the rescaled position and (mi, m2) are parameters chosen 
to give the best spread of lattice points. There should be enough points in the central region 
to accurately follow the rapidly changing fields, while the total number of points should be 
minimized so that the program takes less time to run. 

Once the rescaling is chosen the lattice points are placed equally spread in s between 
s = — 1 and s = 1. This means that the end points of the lattice are at spatial infinity. 

The number of lattice points used was typically between 100 and 300, and the number 
of iterations needed to reach convergence was between 20 and 1000. 

Spatial derivatives are calculated three ways, with forward, backward and symmetric 
derivatives. The energy function E{fi) has contributions from all three forms of the deriva- 
tive, added with equal weighting. The advantage of this approach is that symmetry between 
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left and right is maintained, while odd and even numbered points are kept in close contact, 
which would not be the case were only symmetric derivatives used. 

The values of the fields on the boundaries are set equal to their vacuum values, and are 
kept fixed. The walls are symmetrical between left and right, because of the Zn symmetry 
in the Lagrangian. This links the field values in one half of the wall with the conjugates of 
the values in the other half, i.e. 

S{s) = ^(-s)e'2"'=/^, Y{s) = Y{-s) . (46) 

However, rather than model half of the wall, it was decided to model both sides, because 
of the difficulty of establishing simple boundary conditions at 2; = 0, consistent with the 
strategy of programming with complex, rather than real, fields. 

The convergence criterion was framed in terms of c = \ dE/dfA\- The algorithm is 
judged to have converged when c reaches ^ 10~^°. 

When far from a solution the size of the update may have to be reduced. It is guaranteed 
that, for a sufficiently small movement, following —dfA will approach a solution. However, 
a full NR update may overshoot, and lead to divergent field values, rather than converging. 
To deal with this, in the initial stages of a run the size of the update is reduced by a factor 
of 10, or even 100. Systematic methods of doing this, such as backtracking along dfA, were 
judged to be unnecessary. 

To check the BPS saturation of a solution two criteria are used. Firstly, the total energy 
of the wall is compared to the BPS energy, which is itself easily calculated. Because of 
discretisation error, the numerically determined energy is usually below the true BPS value. 

The second method of determining whether solutions are BPS is to check that the wall 
follows a straight line in W^-space. The straightness of the line is quantified by the ratio of 
the maximum deviation from a straight line linking the ends, divided by the length of the 
line. If this quantity can be shown to decrease towards zero as the number of lattice points 
is increased, the configuration is taken to be BPS. 



B Detailed analysis of the BPS constraint equations 

In this Appendix we supply details of the analysis behind Section ^, where we described the 
values of N and k admitting BPS domain walls. 



In Appendix |B.1| we study the limit A — > 0, where our models tend to the Landau- 
Ginzburg models LGl and LG2 described in Section ^ In this limit we are able to prove 
exact results giving the values of (A^, k) for which BPS walls exist and, conversely, those for 
which BPS walls do not exist. Our argument uses continuity in A to make the statement 
that BPS walls cannot exist in a particular model at A > if its 2; = constraint curve, 
Eqs. ( |3T| , |32D , reaches the LG limit at a point where we know there is no BPS domain wall. 

To this end, we must analyse the z = curves, to enumerate all possible BPS domain 



wall solutions, and to check that they have an LG limit. This is done in Appendix [B.2 
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B.l Analysis in the Landau— Ginzburg limit 

In the Landau-Ginzburg limit, A ^ 0, the extra field Y stays at its vacuum value Y^, = A^=pl. 
We require that solutions of the constraint equations Eqs p^j30| ) exist for all /3 in the range 

< f3 < 27rk/N, and asymptote to i? = 1 at /5 = and (3 = 27Tk/N. Any solution to 
the constraint equation at z = which does not connect the vacua is a 'fake' solution and 
cannot represent a BPS domain wall. By continuity in A, solutions to the full constraint 
equation for our interpolating models which tend to these fake solutions as A — > are also 
inadmissable as BPS walls. 

Firstly, consider the Model I constraint equation Eq. (^), and define a new angle S = 
P — Tfk/N, so that the constraint function becomes 

Clgi{R,S) = R'-^cos[{N-1)S]{-1)'' + {N-1)RcosS-Ncos(^^^ . (47) 

Where cos[(A^ — 1)5] (—1)^ > 0, the constraint function has a single minimum as a function 

01 R at R = -Rmin- One can show that as one varies 5, the maximum value of this minimum 
occurs at 5max.m = Tcm/N, with —k < m < k, and m even (odd) when k is even (odd). At 
these saddle points 

cos iy-Jj-J ~ ( jV ) ) • 

When this is greater than zero there can be no solutions to Clgi{.Ri ^) = 0. For k > 1, there 
is always a value of m for which CiGiiRmin, 5max,m) > 0, and so the constraint curve cannot 
move continuously between 5 = and 6 = nk/N. 

Hence we conclude there are no BPS domain walls in the Landau-Ginzburg limit (A — > 0) 
of Model I for any k > 1, and only one for k = 1 

For the Landau-Ginzburg limit of Model II, we can rewrite the constraint function as 

Clg2{R, 5) = R^^\os[{N + l)5]{-lf-iN + l)Rcos5 + N cos (^^^ . (49) 

Where cos[(A^ + 1)5] (—1)'^ > the constraint function has a single minimum at i? = -Rmin- 
One can again show that the maximum value of this minimum occurs at -Rmin = 1, 5max,m = 
7rm/N, with —k < m < k, and m even (odd) when k is even (odd). At these saddle points 

CLG2{Rmm, 5max) = N (^COS (^^^ - COS (-^)^ ■ (50) 

In this case, the value of the constraint function at the saddles is always less than zero 
(unless m = k), and therefore there are two solutions to the equation Clg2 = when 
cos(7rA;/A^) > 0, and only one if cos(7rfc/A^) < 0. Where cos[(A^ + 1)5] (-1)^= < 0, the 
constraint function becomes a monotonically decreasing function of R. For these values of 
5 there can only be one solution to the constraint equation for coslnk/N) > 0, and none 
for cos{nk/N) < 0. Hence only one of the solutions to the constraint equation at 5 = 
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(which one can straightforwardly see is the one with smaller R) can connect two vacua for 
cos{7ik/N) > 0, and none for cos{7ik/N) < 0. 

In fact, we can also eliminate the solution at = for cos{Tck/N) = 0, as it cannot 
interpolate between vacua with R = 1. 

Hence we conclude there are no BPS domain walls in the Landau-Ginzburg limit (A — > 0) 
of Model II for any k > N/2, and only one for k < N/2 



B.2 Analysis of constraint curves at z = 

In this Appendix we show that the constraint curves aX z = fall into a small number of 
classes: four if N is odd, and six if N is even. These classes are distinguished by having the 
same number of intercepts with the lines representing the LG and VY limits, and by having 
the same number of points at which they are either horizontal or vertical. This classification 
helps us to be certain that we have found all possible types of BPS wall solution in Section |^. 

Firstly we note that, as k is an integer, cos(fc7r) can only have the values 1 or —1 and, 
if N is odd, cos{kn/N) can be either positive or negative. There are therefore four distinct 
cases for which we can expect the constraint curve to be qualitatively different. If is even, 
cos{kTT/N) can be zero, and there are therefore six classes. 

We have already eliminated cos{TTk/N) < from consideration in the Landau-Ginzburg 
limits of both models, and thus we need concentrate only on the cases cosvrfc = ±1, with 
cos{7rk/N) > For illustration purposes it is therefore sufficient to consider, as we have 
done, N = 5 with k = 1,2 only. 

We look for intercepts with four particular curves, namely 

i?o = 1, (51) 

Yo = 0, (52) 

Yo = NtI, (53) 

Yo = NtI- NlnRo , (54) 

with the negative sign taken for Model I and the positive sign for Model II. Note that the 
third of these is the LG limit, and the fourth is the VY limit in the case of even k. Along 
these curves the constraint equations are sufficiently simple that the number of solutions can 
easily be deduced. 

The other information we use is the slope of the constraint curves 



dYp 

ORn 



(iV-l)i.-e-o..-.(_i)._^„ Model I 



i?o(^o'^e-^o+iv-i(_i)fc _ 1) 
(Ar + l)i?^e^o-A^-i(_l)fc_yp 



Ro{Ro e 



NpYo-N-l( 



(55) 

Model II 



By noting the sign of dYo/dRo we can trace the shape of the curve. Note that for even k 
the constraint curve is vertical at Yq = N ^1 — NIuRq, which is precisely at the VY limit. 
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This is illustrated in Figs 5b and 5d, where it can be verified that the constraint curves are 
vertical at Rq = Rqvy- 

Now we shall consider the two cases {k even or odd) in detail, in Model I. 

When k is even (see Fig 5b), Eqs. (^,^B|) can not be satisfied, so the constraint curve 
never crosses the lines i?o = 1, = ^ =F 1- However, it does cross Eq. (0) once at -Rq > 
and Eq. (|5^ ) twice, once on either side of i?o = 1- 

Let us trace the curve from its Ro < 1 intercept with Eq. (0). For the part of the curve 
with Yq > N — 1 — NIiiRq, the slope is negative, and so the curve must move left and 
up. The part with Yq < N — 1 — NIiiRq has positive slope but, because it cannot cross 
Yq = N — 1, it must reach a turning point, and also move left and up. This accounts for the 
upper branch of the curve in Fig. 5b. 

Starting now from the Rq > 1 intercept with Eq. (|5^), the lower part of the curve will 
follow Rq oo,Yq — oo. The upper part will cross Eq. and then BYq/ORq = and 
so turn down again. It cannot cross Eq. (|52|) again, and so Iq ^ from above, as Rq oo. 

When k is odd, (Fig 5a), we can by inspection of Eq. ( |5^ see that Ci{Ro,Yq) = only 
when Yq > 0. Furthermore Eqs. ( pI| , |53D can each be satisfied once, while Eq. (|5^) is never 
true: hence the curve is never vertical, and in fact the slope is always negative. This means 
that lo from above as Rq —* oo, while as Rq —>■ 0, Yq oo. 

One can go through a similar analysis to satisfy oneself that the curves for Model II are 
similar to Figs. 5c, 5d. The principle differences between the models occur for even k and 
can be seen by comparing Fig. 5b to Fig. 5d: Model I has no intercepts with Rq = 1, while 
Model II has two. There are also no intercepts with Yq = N — 1 in Model I, but one can easily 
check that there are two with the analogous line Yq = N + 1 in Model II. The intercept with 
J?o > 1 does not however correspond to a domain wall in the Landau-Ginzburg (A —>■ 0) 
limit, as was shown in Appendix p.l| . 
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